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Abstract 

We study a pair of commuting difference operators arising from the elliptic solution of the 
dynamical Yang-Baxter equation of type Ci. The operators act on the space of meromorphic 
functions on the weight space of sp(4, C). We show that these operators can be identified with 
the system by van Diejen and by Komori-Hikami with special parameters. It turns out that 
our case can be related to the difference Lame operator (two-body Ruijsenaars operator) and 
thereby we diagonalize the system on the finite dimensional space spanned by the level one 
characters of the C^-affine Lie algebra. 

1 Introduction 

The Ruijsenaars system of difference operators [R], are difference analogue of the Calogero-Moser 
integrable system of differential operators. The operators of the system is denned in terms of 
elliptic function, and in the trigonometric limit, they degenerate to the Macdonald ^-difference 
operators [M]. The Ruijsenaars system has been studied extensively. Especially, Hasegawa shows 
that this system can be obtained as transfer matrices associated to the Sklyanin algebra [H3] and 
Felder-Varchenko reconstructed them as transfer matrices associated to the dynamical i?-matrices 
[FV2]. These two approaches are related by the vertex-IRF correspondence [Bax] [JMOl]. 

Extending these works, in [HIK] we construct a pair of commuting difference operators acting 
on the space of functions on the C2 type weight space. The method therein is based on the elliptic 
solution of the dynamical Yang-Baxter equation of type Ci (or Boltzmann weights of the face 
model [JM02]). We have also shown that the space spanned by the level one characters of the 
affine Lie algebra sp(4,C) is invariant under the action of the difference operators. 

On the other hand, a generalization of the Ruijsenaars system to BC n case is studied by 
van Diejen [vDl] and Komori-Hikami [KH1] [KH2]. First, van Diejen constructed two elliptic 
commuting operators, one is of the first order and the other is of the n-th order. Therefore he 
obtained an elliptic extension of difference Calogero-Moser system of type BC2 [vD2]. Extending 
this work by van Diejen, Komori and Hikami obtained a general family of n commuting difference 
operators with elliptic function coefficients. Besides the step parameter of difference operator and 
the modulus of elliptic functions, the family contains nine arbitrary parameters. Their construction 
uses Shibukawa-Ueno's elliptic i?-operator [SU] together with the elliptic X-operators [KH3] [KH4], 
the elliptic solution to the reflection equation. It can be regarded as an elliptic generalization 
of Dunkl type operator approach to Macdonald systems, which have been extensively used by 
Cherednik [Ch] (see [N] for BC n case). 
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This paper has two goals. One is to establish the relationship between our system of difference 
operators and the van Diejen-Komori-Hikami system. The other is to diagonalize our difference 
operators on the finite dimensional space spanned by theta functions. The first goal is attained in 
section 2 and second is in section 3. 

In section 2 we review the construction of the elliptic difference system of type C 2 and give a 
new form of our operators. After this, we will establish an identity consisting of theta functions 
(Lemma 1), and explain how our system can be identified with from van Diejen-Komori-Hikami 
system with special choice of parameters (Theorem 2). That is, our approach to the difference 
operators as transfer matrices, based on the knowledge of the Boltzmann weights, reproduces 
a special case among the family of commuting operators obtained by Dunkl-type approach. It 
should be also mentioned that those two approaches to the system are not yet related, although 
the resulting commuting operators have the relationship as above. 

In section 3, we introduce the finite dimensional space of theta functions invariant under the ac- 
tion of Weyl group and its basis after Kac-Peterson [KP] . Our aim is to diagonalize our operators on 
this space (Theorem 4) . This is an elliptic analogue of the eigenvalue problem of Macdonald opera- 
tors on the space of symmetric polynomials. Their eigenfunctions, called Macdonald-Koornwindcr 
polynomials, are much investigated in g-orthogonal polynomial theory [K] [N]. 



2 The difference operators of type C2 

2.1 Construction of the difference operators of type C2 

Let 2 be the Lie algebra sp(4, C), t) its Cartan subalgebra and f)* the dual space of 1). We realize 
the root system R for (g, f)) as R := {±(ei ±£ 2 ), ±2£i, ±2e 2 } C f)*. A normalized Killing form ( , ) 
is given by 



(sj,e k ) - -Sjk, 



(2.1) 



and the square length of the long roots ±2ei is two. We will identify the space f) and its dual [)* 
via the form ( , ). The fundamental weights are given by Ai = £1, A 2 = £1 + £ 2 . Let Vd be the set 
of weights for the fundamental representation L(Ad). We have 



Vx = {±£!,±£ 2 }, V 2 = {±(£l±£ 2 ),0}. 



Let d, d' be 1 or 2. The type Boltzmann weights of type (d, d') are given as follows. Fix a 
complex parameter fteC. For any pair A, /j,,v,k (z f)* of weights, the Boltzmann weight 

W dd> 



(2.2) 
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This equation is also known as the dynamical Yang-Baxter equation. Here we give the explicit 
formula for W n and see [HIK] for the other type W dd > {(d,d') = (1,2), (2,1), (2,2)) which are 
obtained by fusion procedure. They are expressed by the Jacobi theta function 9\{u) — 6\{u\t) 
with elliptic modulus r in the upper half plane fj + (See Appendix B for the definition of 9\(u)). 
For p, q, r, s G V such that p + q = r + s, we will write 



P 



A 



A + 2hp 



s^\q= Wu ^ x + 2hg x + 2h(p + q) 

The explicit formula for Wn is given as follows: 
?i(c-u)0i(u + ft) 



V 

p\u\p 
P 



P 

p\u\q ■ 

q 

Q 

p\u\p ■ 

q 

q 

p\u\-q- 

-P 

P 

p\u]-p 

-P 



0i(c)0i(ft) ' 

6>i(c - u) 9i{\ p - q - u) 
0i (c) 0i(A p _ 9 ) 

9 l {c^u)9 l {u)9 l {\ p ^ q + K) 



(p + ±q), 
(p ¥= ±q), 



9i(c) 0i (ft) 0i (V*) 

01 (u) 01 ( A p+g + ft + C - U) 01 (2 A p + 2ft) Ur^±p gl ( A P+r + ft) 

0i(c)0 1 (A p+9 + fi) 0i(2A ? ) n^±,«i(Vr) 

0i(c-u) 0i(2A p + ft-it) 
0i (c) 0i(2A p + ft) 

?i (u) 0i (2 A p + ft + c - m) 0i (2 A p + 2ft) -,-r 0i(A p+g + ft) 



(2.4) 
(2.5) 
(2.6) 

(P^«), (2-7) 
(2.8) 



0i (c) 0i(2A p + fi) 



i(2Ap) 



n 



l(A p + 9 ) 



Here the crossing parameter c is fixed to be c := —3ft. 

We define the difference operators M d (u) (u G C,d = 1,2) acting on the functions on \) by 
means of the Boltzmann weights of type (1, 2) and (2, 2). 



(M d (u)f)(X) := ]T W d2 
pev d 

Here the shift operator T$ is defined as 



A A + 2hp 
A A + 2hp 



u) TU{\). 



Tlf(X) :=f(\ + 2hp). 



For A 6 I)* and p G P d (d = 1, 2), we put 



(A,p). 



Note that if we denote Aj = (A,£j) (i = 1,2) and /(A) = /(Ai, A2), then 

T* 2e J(\ u \ 2 ) = /(Ai ± ft, A 2 ), T^ 2£2 /(Ai, A 2 ) = /(Ai, A 2 ± ft). 
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Theorem 1 [HIK] 

(i) For each u, v 6 C, we have M d (u)M d ,(v) = M d ,(v)M d (u) {d,d! = 1,2). 

(ii) The explicit form of M d (u) are as follows : 

M l{ u) = F(u) £ J] °%T L ^r T l, (2-9) 



p=±ei 
\q=±E2 



(2.10) 



Here U(X P , X q ) is given by : 



u(x p ,x q ) = -± 



(2ft) 0i(2A p + 2ft) X (2A 9 + 2ft) 0i(A p+g - 5ft) 0i(A p+9 + 2?l ) 



0i(6ft) 0i(2Ap)0i(2A 9 ) 0i(Ap +9 )0i(Ap +9 + ft) 

and -F(u), G{u),H(u) are the following functions depend only on u and ft : 



F(u) := — 

ft, („ - 

G{u) 



0i(-3ft) 2 0i(ft) 2 ~ ' 
0i(u - ft) 9 l (u) 2 6>i(u + ft) 0i (u + 2ft) 0i (u + 3ft) 2 0i (u + 4ft) 



and -ff(w) := 



0i(-3ft) 4 0i(ft) 4 

0i(u + 6ft) 0i (u- 3ft) 0i (2ft) 
0i (u) 0i(u + 3ft) 0i (6ft) 



The following Lemma is the key for the identification with van Diejen's system as well as 
for the diagonalization of our difference operators. The author is grateful to van Diejen for the 
information. 



Lemma 1 We have 



Ett(\ \ \ ^i(A P+g - ft) 0i(Ap +g + 2ft) 

U(Xp,X q ) - ^ 0l( Ap +9 )0 l( Ap +9 + ft) ~ K > (2 - 11} 



p=±£i p=±ei 

g=±S2 9=±£2 



wftere K is a constant given by 



(8ft) 0! (ft) 0! (5ft) 0i (2ft) 0i (6ft) 0! (3ft) 0i (4ft) 0i (ft) 



0i (6ft) 0i (5ft) 0i (4ft) 0i (3ft) 0i (5ft) 0i (4ft) 0i (3ft) 0i (2ft) ■ 

Proof. Let /(A p ) be the left-hand side of (2.11), regarded as a function of A p (p € F). It is doubly 
periodic function of the periods l,r. Let us show that it is entire. The apparent poles of f(X p ) are 
located at 

Ap = X q , Ap = X q - ft (p,q e I, p + q ^ 0), Ap = (p e F). 

Note that /(A p ) is TT-invariant, then the points A p = X q and A p = are regular. Also, the residue 
of /(Ap) at Ap = — X q — ft is 

0i(2ft)0i(-2A 9 )0 1 (2A (? + 2ft)0 1 (-6ft)0i(ft) 0i(-2ft) 0i(ft) _ () 



0i (6ft) 0i(-2A g - 2ft)0i(2A 9 ) 0i (-ft) 0i(-ft) 
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Now we have proved that f(X p ) is independent of A p , then we consider g(X g ) = /(— X q — 2ft) as 
a function of X q (q ^ p e I): 

9i (2ft) ( 0i (2A g + 2ft) 6>i (2A g - 2ft) 6» 1 (2A <? + 7ft) 

9\Xq) - 



0i (6ft) V 0i (2A 9 + 4ft) 0i (2 A 9 + 2ft) 0i (2A 9 + ft) 
i(2A g + 6ft) 0i(2A 9 + 2ft) 0i(2A 9 - 3ft) 



0i(2A g ) 0i(2A, + 2ft) 0i(2A g + 3ft) 
9i(2A 9 + 6ft) 0i(2A 9 - 2ft) x (-3ft) 0i(4ft) 
0i(2A 9 + 4ft) 0i(2A g ) 0i (2ft) 0i (3ft) 
0i(-2Ag-3ft)0i(-2Ag) 0i(2A 9 + ft)0i(2A (? + 4ft) 0i(ft)0i(4ft) 



0i(-2A g -2ft)0i(-2A g -ft) 0i(2Ag + 2ft)0i(2A 9 + 3ft) 0i(2ft) 0i(3ft) ' 

By the same argument we can show that g{X q ) is independent of A g . Therefore we get K by putting 
X q = ft in g(X q ) and the proof completes. □ 

2.2 Identification with van Diejen's system 

We define the difference operators Md to be the components of Md{u) independent of u: 

~ = E n ^_A T ^ (2 . 12) 



~ ^ ( Oi{X p+q -K) h h Oi(X p+q - h)6 1 (X p+q + 2h) \ 

p=±ei 

9=±£2 



0!(A p + g+ft) ^ 0i(A p + g)0i(A p + g+ft) 



More general commuting difference operators Hi,H 2 are obtained by van Diejen_and_ later by 
Komori-Hikami in a different way. In this subsection we identify our operators Mi,M 2 as van 
Diejen's system of difference operators with special values of parameters. The operators Hi,H 2 
depend on nine complex parameters fi, fi r , \j! r (r = 0, 1, 2, 3) satisfying the condition 

^( Mr + M ;) = (2.14) 

r 

and are defined by 

Hi = ^ w(exi)v(exi + x 2 )v(ex 1 — x 2 )T] 1 
ti 

w(ex 2 )v(ex 2 + xi)v(ex 2 - X!)T^ 2 + [7{i,2},i, 



e=±l 



e=±l 



H 2 = w{exi)w{e'x 2 )v{exi+e'x 2 )v{exi+e'x2+ r y)T] 1 T] l2 

e,s'=±l 

+ U{ 2 -^! ^2 w(exi)v(exi + x 2 )v(ex! - x 2 )T] 1 

e=±l 

+ E w ( £a; 2)w(£a;2 + xi)v(ex 2 - xi)T] 2 + f/{i,2},2- 



e=±l 
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Here (i = 1, 2) stand for the shift operators 

T£ 1 f(x 1 ,x 2 ) = f(xi ±7,x 2 ), Tl 2 f( Xl ,x 2 ) = f(x 1 ,x 2 ±'y) 

and 

where cr(z) = <7o(;z) denotes the sigma function with two quasi periods u)\, uj 2 and a r (z) (r = 1, 2, 3) 
associated function obtained by shift of argument over the half periods (See Appendix B for more 
detail). The functions J7{j},i, ^{i,2},j (j = 1,2) are defined as follows : 

U {J}A = - w ( x j) - w (- x j) U = 1. 2), 

\T- TT Oy(M - 7/2 + Xj) oy(/i - 7/2 - Xj ) 
{1 ' 2},1 ~oks ^(-7/2 + ^)a r (- 7 /2-^) ' 



where 



2 

* = gfrM/i - 7) II ^(/V W -7/2K(/C( s) ), 



0<s<3 

with 7r r denoting the permutation 7r = id, ii\ = (01)(23), ir 2 = (02)(13), 7r 3 = (03)(12). 

U{i.2},2— ^ w(exi)w(e'x 2 )v(exi + e'x 2 )v(—exi — e'x 2 — j) (2-16) 

£,£'£{1,-1} 

We mention that the Komori-Hikami system in [KH2] is of more complicated form and has nine 
arbitrary parameters, that is, they removed the condition (2.14). 

In Hi,Ti. 2 , we specialize parameters /x, (r = 0,1,2,3) as /1 = —7, \i r — \j! t — 0. Then 
w(z) — 1 and [7(1,2}, 1 = 0. Let us denote these specialized operators by Hi,H 2 . Because of these 
simplifications^ we^jmmediately obtain the following from Lemma 1, giving the identification of 
our system {Mi, M 2 } and van Diejen's {Hi, H 2 }- 

Theorem 2 For a function f(Xi,X 2 ) on f), we set (p(f)(x\,x 2 ) by 

its, s mixl + xl) x 1 x 2 

<p(f)(xi,x 2 ) :=exp i ^/( ), 

Wl Zul\ llO\ 

and let 7 = 2uih, we have 

ip Mi ip- 1 = e 2 "iT 2 /^i 7^ 

if M 2 ^ = e 2 "iT 2 /^i {fc 2 + 2-Hx) . 

Proof . Use the connection between the theta function and sigma function (B.8) in Appendix B 
and (2.11) to compare (2.10) and (2.16). □ 
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3 Diagonalization of the system 
3.1 The space of theta functions 

Let Q and Q v be the root and coroot lattice, P and P v the weight and coweight lattice respectively. 
Under the identification i) — i)* via the form ( , ), they are given by 



and 



P = ^ Zej, Q v = Z2ej 

3=1,2 o = l,2 



P y = Q = Z2 £l + Z2e 2 + Z(ei + e 2 ). 



(3.1) 



For (3 G 1)*, we introduce the following operators T T p, Tp acting on the functions on ()*: 
(7>/)(A) :=/(A + /3), 



(T T/3 /)(A) :=exp 



2m ( ( A ,/3) + Mr 



/(A + r/3) 



We define the space of theta functions (of level 1) by 

T/ii := {/ is holomorphic on fj* | T Ta f = T a f = f (Va G Q v )} . 
For each fj, G P and fixed r G we define the classical theta function 6 M (A) of A <E f)* by 



2™ ( ( 7 ,A) + %^r 



It is known that 



e ^( A ) : = exp 



{e M (A) |/i = 0, ei, e 2 , £i + £2 modQ v } 



gives a basis for T/ii over C [KP]. 

Let W C GL(f)*) denote the Weyl group for (g, f)), and consider the M^-invariants in T/ii: 

Thf := {/ € T/n | /M) = /(A) (Vw G W)} . 
Theorem 3 [HIK] The operators M x , M 2 preserves Thf . 

For [i G P , we define W" M := {w G W | w/i = /i} and introduce the following symmetric sum of 
theta functions, 



5 a»( A ) : - jTFy H e «.(p)( A )- 



Then 



{S P (A) a/ = 0, Ai (= ei), A 2 (= £l +£ 2 )modQ v } 

forms a basis for Tlif over C. 

It is known that ThY is also spanned by the level 1 characters of the affine Lie algebra Sp(4, C). 
Note that 0_ M (A) = M (A) and £l+£2 (A) = £l _ £2 (A). So that we have 



5 (A) = e (A), 5 Al (A) = 2(e ei (A) + e2 (A)), S A2 (A) = 40 £l+£2 (A). 
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3.2 Diagonalization of M d 

In this subsection, we diagonalize the operators M d on the space ThY ■ We set 

/i(A) :=6 £l (A)+e £2 (A), / 2 (A) :=e (A) + e £l+£2 (A) and 
/ 3 (A) :=e (A)-9 £l+£2 (A). 

They are linearly independent in the space ThY ■ 

Theorem 4 The functions fi(\) (i = 1,2,3) are common eigenf unctions of Md : 

M d fi(\) = E d>i fi(X) (d =l,2,i = 1, 2, 3). 
The eigenvalues are given by 

M(2fi)fl m (0) \ 2 

and E 2 .i = 2£i,i, where the Jacobi theta functions 8i(z) = 6i{z\r) (i = 2,3,4) are defined as in 
Appendix B. 

We willjjrove this theorem by using the following three lemmas. First, we show that the 
operators Md split into two Ai-type components. 

Lemma 2 Let us denote A± := (A, £i ±e 2 ) and define 



_ 3i(\±-h) n fli(-A±-ft) h 
± ' 0i(A±) £1±£2 0i(-A±) -( £ i±^)- 



Then we have 

Mi = M 2 = H 2 + + H 2 _. (3.2) 

Proof . To prove the first identity, we note that 

0i(A+-ft) h 1 (A_-ft) » 



(A+) ^ £l+£2 0i(A_) ei ~ £2 



7 H A + 

fli(A+-/Q gi((A + ft(ei+ g2 ))--fi) s ft 
0i(A+) fli((A + ft(ei+e 2 ))-) £1+£2 £1 ~ £2 



t(A+-ft) ^(A_-ft)^ 



0i(A+) fli(A_) 

Here we used the identity (ei + e 2 , £i — e 2 ) = 0. The second identity follows from, for instance, 

gi(A+-j) » 0i(A + -fi) fi 

0i(A+) £l+£2 0i(A+) £l+£2 
= 0i(A + -fi) g 1 ((A + fi(g 1 + £2 )) + -fi) s ft 
(9i(A+) (9i((A + ft(ei+e 2 ))+) £l+£2 £l+£2 
0i(A + -fi) g 1 (A++fi-fi) ft 
(9i(A+) 0i(A+ + /i) 2 ( £l+£2 ) 
_ #i(A+ - ft)„R „ft 



?i(A+ + /i) ^ 



2ei J -2e 2 - 
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Here we used the identity (si + £2, £1 + £2) = 1- □ 

Second, we consider the eigenvalue problem for the Ai-type difference operator (difference Lame 
or two-body Ruijscnaars operator) 



(z-ih) 6i(z + eh) 



f(z + h)+ ° iy :," l> f{z -h)= Ef(z). (3.3) 



0i(z) JK ' Otiz) 
Lemma 3 For the special coupling constant £ = 1, the functions 

9i(z) (* = 2,3,4) 
are solutions of the equation (3.3) with eigenvalues 

0! (2ft) 0,(0) 

E - Ei - mm (i - 2 ' 3 ' 4) - 



Proof. We note that the functions 02 (z), 03 (.2)7 and 6i{z) can be rewritten in e cz 9i(z + 1) up to a 
constant, where 

(t,c) = (^,0), (^.Tri) and (I Tri) (3.4) 

respectively (See the formula (B.3) in Appendix B). Let (t,c) be one of these, and we denote by 
g(z) the function which obtained by the action of the Ai-type difference operator (3.3) with £ = 1 
to e cz 9 x {z + t): 

_ e^z-h) 0i(z + t + h) C(z+R) e 1 (z + h)e 1 (z + t-h) (z _ h) 

9{) - 9i(z) + 

This is holomorphic and doubly quasi-periodic function: 

g(z + 1) = -e c g(z), g(z + r) = -e™ - 2m{z+t)+CT g(z) . 

Moreover, g(z) — at z = —t. Therefore, g(z) is equal to e cz 9i(z + i) up to a constant, which is 
the value of 

9 1 (z-h)9 1 (z + t + h) h 9 1 (z + h)9 1 (z + t-h) h 

Bx{z)Bx(z + t) 0i(z)0i{z + t) 1 ' ' 

at any chosen point. If we choose z = h, then the first term in (3.5) vanishes and we have 

0i (2ft) 0i (t) _ ch 0i(2ft)0 i (O) 



0i (ft) 0i (ft + t) 0i (ft) 9, (ft)' 

where i = 2, 3 and 4 corresponding to the values of (t, c) in (3.4), as an eigenvalue. □ 

Remark. This can be regarded as a special case of Felder-Varchenko's study [FV1]. They 
expressed the solutions of (3.3) in terms of the algebraic Bethe Ansatz method, which is originally 
developed and applied to the spin chain model. In fact, the operator in the left hand side of (3.3) 
can be regarded as the transfer matrix of the simplest spin chain, that is, it consists of only one 
site of freedom with spin £ = 1. In this case, the Bethe Ansatz equation 

0l(*-ft) _ 2Hc 



L (i + ft) 



e 2Rc (3.6) 
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is exactly the same as the condition that the function (3.5) dose not have pole at z = —t. 
Because of Lemma 3, the product of the theta functions 

^(A+)^(A_) (i,j = 2,3,4) 

are simultaneous eigenf unctions of the operators H+,H?_ and H + H- with eigenvalues 

e^hfem 2 e^nf ej(o) 2 e^hf ^(o)^(o) 
e^h) 2 6i(h) 2 ' e^hfe^h) 2 an e 1 {h) 2 9 l {h)0 1 (h) 

respectively. Finally, we shall establish the relationship of these Bethe Ansatz solutions and the 
bases of Thf. 

Lemma 4 The functions f%{\) £ ThY are expressed in terms of the Jacobi theta functions as 
follows: 

/i(A) = 2 (\ + )e 2 (\_), / 2 (A) = 3 (A+)0 3 (A_), / 3 (A) = 4 (A+)0 4 (A_). 

Proof. Because of the definitions of coroot lattice Q y (3.1) and Killing form (2.1), each basis of 
Thi is expressed as 

6 (A) = 3 (2Ai|2t)0 3 (2A 2 |2t), 
ei (A) = 3 (2Ai|2t)0 2 (2A 2 |2t), 
6 £2 (A) -0 2 (2A 1 |2t)0 3 (2A 2 |2t), 
9 £1+£2 (A) =0 2 (2A 1 |2t)0 2 (2A 2 |2t). 

Here A^ = A £i (i — 1,2). Therefore we can prove this lemma by using the identities of theta 
functions (addition theorems) (B.4), (B.5), (B.6), (B.7) in appendix. □ 

We note that the anti-symmetric function 6 £l (A) — £2 (A) = #i(A + )#i(A_) is also the eigen- 
function with eigenvalue zero. 
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Appendix A: Differential limit 

Let us clarify the connection between our system of difference operators and a quantization of the 
Inozemtsev Hamiltonian [IM] [I] . By expanding in ft one infers that 

Mi = 4 + M!, 2 ft 2 + M M ft 4 + 0(ft 5 ), 
M 2 = 8 + M 2 , 2 ft 2 + M 2j4 ft 4 + 0(ft 5 ). 
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If we abbreviate a function /(A £l ± E2 ) as /(±), di — (i — 1,2), and 0[(z) = -^6\{z) etc. We 
have 



Mi, 2 + ^ 

-2 (£(+) + -2(^( + ,-^(-)U 



'/)// a// 



M 2i2 =2Mi >2 , 



and 

M 2i4 - 2M M = d 2 <9f 



- 2 (^ +, -|<->) af *- 2 (l (+,+ l ( - ) ) aA2 

l) 2,+)+ (l) 2, - ) )-(f (+)+ 4 (+ »l < -» + f<- , )l 3f 



£) 2 M£) 2( - , )-(? <+, -4< + >£<-> + £<- 



+J|('£V(+>-(T) 2 (->W 



+ { 2 (f< + ) + f(-)) +2 (f( + )|(-) + |( + .f(-) 

&)'<♦>♦(*)'<->)}* 
^<+>-^h)- 2 (£(+$(-) -£<+)£(-> 

+ HV (+)+ V<-»)- 4 (^< + » + t ( - ) 

(/)///)/ 2 2 \ /J// /J// 

We set A = (9i(+)6>i(— ), then 

A-.M 2 , 2 .A = ^^I + 4((f-f)( +)+ (f-f)(-) 

= 9 1 2 + a 2 2 + 4((log0 1 )"(+) + (logei)"H), (A.l) 
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A" 1 • (M 2;4 - 2Mi. 4 ) • A 
dldl 



" Ql 2 \ / /)// /)/ 2 



+2 



0i ; \0i ej 

,,9'" JlQ'l J'l 3 \ , x /^l" ^i 3 \ , A, 

- ((-?-- 3^ + 2 ir ) (+) + 2 ^J (-)J 92 

lllnl Qlllnl 




9? 

(-) I -*£<+)£(-) 

1W1 H)+8(^(+)^(-) + g-(+)^(-) 

/ \ 2 / qi \2 

= 

+4{(log0i)"(+)-(log0 1 )"(-)}0 1 fc 

+2{(logfl 1 )"'(+) + (log^n-)}^ + 2{(log0 1 )'"(+) - (log^)'"(-)}a 2 
+2{(logfl 1 )( 4 )(+) + (logfl 1 )( 4 )(-)} 
+4{(log0 1 )"(+)-(log0 1 )"(-)} 2 
= {ftft + 2 ((k>g0i)"(+) - (log0i)"(-))} 2 

The complete integrable Hamiltonian of type BC n is introduced by Olshanetsky-Perelomov 
[OP], and later generated by Inozemtsev-Meshcheryakov [IM] [I]. In the rank two case, the Hamil- 
tonian is 

H = - ^{dl + dl) + g(g - 1) (p(zi + x 2 ) + p{x x - x 2 )) 



X] 9r(9r - l)(p(w r + p(w r +X 2 )), 



0<r<3 



where p(x) denotes the Weierstrass p-function with two periods 2wi and 2w 2 , and wo = 0, w 3 = 
— wi — Ld2- By the connection between theta function and p function (B.9) in Appendix B, our 
differential limit (A.l) is identified with this Hamiltonian for the special coupling constants g(g — 
1) = 2, and g r (g r - 1) = (0 < r < 3). 

Appendix B: Theta function 

We establish notations and identities on the theta functions [WW] . The Jacobi theta functions are 
defined for t £ f) + as follows: 

6Mt) = E CX P [ 2 ™ (( z +\)(k + \) + \(k+ \fr) 
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02(z|t) = y^exp 

fcez 

#3(z|t) = J^exp 



•_>..■ .■ |<fc + i) + i(fc + i) 2 r 



fcez 



27TI ( ZK + — T 



r ( i fc 2 \ 

04(z\t) =22 e *P 2« ( (z + -)fc+yrj 



Note that 9\{z) is odd and the other three are even. These functions has quasi-periodicity: 

e 1 (z + m\T) = (-l) m e 1 (z\r), (B.l) 
9 Y (z + mr|r) = (-^"e-^™ 2 ^ 2 ™^ (z|r), (B.2) 
(m £ Z), while other three can be expressed by 9\{z) 

e 1 {z + l -\r) = e 2 (z\r), 

6 1 (z + ^\T)=ie-*< z+ i)6 i (z\T), (B.3) 

e 1 (z+ l - + T -\r)=e-<^)e,{z\r). 
We use these identities in the computations in Lemma 4. 

# 4 Mt)0 4 (2/|t) = 9 3 {x + y\2T)e 3 (x - y\2r) - 9 2 (x + y\2 T )9 2 (x - y\2r), (B.4) 

Ux\T)9 3 (y\T) = 9 3 (x + y\2r)9 3 {x - y\2r) + 9 2 (x + y\2r)9 2 (x - y\2r), (B.5) 

92(x\r)9 2 (y\T) = 9 3 {x + y\2r)9 2 (x - y\2r) + 9 2 (x + y\2r)9 3 {x - y\2r), (B.6) 

9 1 (x\T)9 1 (y\r) =9 3 (x + y\2r)9 2 (x - y\2r) - 9 2 (x + V \2t)9 3 {x - y\2r). (B.7) 

The sigma function a(z) is an entire, odd, and quasi-periodic function with two primitive 
quasi-periods 2o>i, 2w 2 . 

a(z + 2nuj 1 + 2mLJ 2 ) = (^-l^+m+nm^n^m^Kz+n^+m^)^ 

with r}i = C(u>i) (i = 1,2), where ((z) = cr'(z)/a(z) denotes the Weierstrass (-function. The 
connection between the Jacobi theta functions and the sigma functions are 

j h z 2 \ 0i(*/2wi) 



a(z) = cxp 



2wi 



9[(0) 



a r (z) 



exp 



2wi 



fl r+ i(z/2a>i) 

«r+l(0) 



(r = 1,2,3). 



Then, for the function v(z) in van Diejen's system (2.15), we have 



v(z) :-- 



a(z + (J,) 



cxp ■ 



7? 1 (2z A i + M 2 )\ 9 1 {{z + t l)/2u 1 ) 



2wi 



0i(*/2wi) 



The connection with p function is 



= -£2 1 °g fT ( z ) = -773 



1 

4^ 



k>g0i(*/2wi) 



(B.8) 



(B.9) 
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